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reduced flutter freguency 
lift due to motion of sth blade per 
unit span, positive up 
lift due to wakes of sth blade per 
unit span, positive up 
riondimensional lift coefficients due 
to bending and torsional motions, 
respectively, In rth mode 
nondimensional moment coefficients 
due to bending and torsional 
motions, respectively. In rth mode 
nondimensional lift and moment co- 
efficients, respectively, due to 
wake i'l rth mode 

moment about the elastic axis due to 
motion of sth blade per unit span, 
positive nose up 

moment of sth blade per u’'1t span 
about the elastic axis due to 
wake, positive nose up 
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Mach number 
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guantities defined in Eg. (6) 
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mass per u>Ht span of sth blade 
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bending deflection of sth blade 
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number of blades in cascade 
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bending deflection of blade in rth 

[P] 
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[I] 
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mass moment o'< inertia of sth blade 
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radius of gyration of sth blade. 
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respectively, of sth blade 
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reduced freguency, ub/V 
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static mass moment of sth blade per 

kar 

acoustic resonance reduced freguency 


unit span about elastic axis, 
positive when center of gravity is 


aft of elastic axis 
time 

freestream velocity relative to the 
blade 

flutter speed 

velocity induced by wakes 

column matrix, defined in Eg. (3) 
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rectangular coordinate axes 
dimensionless static unbalance of 
sth blade (-S» /mjb) 
column matrix, defined In Eq. (3) 
ampTitude of torsional motion of sth 
blade, positive clockwise 
torsional amplitude of each blade tf 
tuned cascade 

amolltude of torsional deflection of 
a blaie In rtn mode of a tuned 
cascace 

Interblade phase angle, 2*rl‘n 
Interblade phase angles at acoustic 
resonance 

nondimenslonal eigenvalue, (un/w)‘ 
nondimenslonal uncoupled bending 
frequency of sth blade 
nondimenslonal uncoupled torsional 
frequency of sth blad*' 
logarithmic decrements uf sth blade 
In bending and torsloii, respec- 
t Ively 

damping ratios of sth blade In bend- 
ing and torsion, respectively 
location of elastic axis measured 
from leading edge, (a ♦ l)/2 
mass ratio of sth blade, i*s ■ mj/tpb^ 
real part of eigenvalue, defined In 
Eq. (8) 

imaginary part of eigenvalue, 
defined in Eq. (8) 
nondimensional flutter frequency 
stagger angle. Fig. 1 
fluid density 
frequency 

re ference frequency 
matrices 


(■) 


-1 


differentiation with time 
inverse of a matrix 
indicate summation over r - 0, 1, 
^' ... ( 1-1 


1. INTRODUCTION 


The prediction of aeroelastic stability and re- 
sponse of bladed-disk assemblies used in aircraft 
turbofan engines has been receiving considerable 
attention in the literature. While in most of the 
research in th>s area all the blades are assumed 
identical, a limited amount of workl“^ considering 
either pure bending or torsional motion of the blades 
has also included small differences between individ- 
ual blades, known as mistuning. The results in these 
references have shown a beneficial effect on blade 
flutter and an adverse effect on forced response. In 
spite of these findings, the lack of complete under- 
standing of the phenomena has prevented its in- 
corporation In current design analysis systems. 
Furthermore, the importance of the effects of Inher- 
ent random mistuning due to manufacturing tolerances 
and of controlled mistuning upon the aeroelastic 
characteristics of bladed-disk assemblies have been 
demonstrated by actual engine experince.^*®*®. 

A research program in turbofan engine aero- 
elasticity is being conducted in the NASA-Lewis 
Research Center. As a part of this general program, 
an effort is in progress to improve the basic under- 
standing of turbofan engine aeroelastic characteris- 
tics including mistuning effects. The effects of 
blade mistuning on coupled bending-torsion flutter 
and aeroelastic response due to wakes in incompressi- 


ble flow have been studied In Ref. 10. The purpose 
of th<i paper Is to continue and to extend the work 
of Ref. 10 Into the subsonic and supersonic flow 
reglme.i. To the best of the author*.' knowledge, the 
Aeroelastic characteristics of mistuned cascades In 
subsonic and supersonic flows using a multl-degiee of 
freedom model for the blades have not been studied In 
the published literature. 

Except for the unsteady aerodynamic models; 
the mathematical model considered Is the same as 
that used In Ref. 10. In the present paper, the 
unsteady aerodynamic loads are calculated by using 
Smith's*^ theory In subsonic flow and Adamezyk and 
Goldstein's^^ theory In supersonic flow with a 
subsonic leading edge. A brief discussion Is pre- 
sented on three regimes of fl 'tter, which are based 
on the decaying or propagating nature of an acoustic 
pressure disturbance In both subsonic ai'd supersonic 
flows. 

II. THEORY 

A. Structural Model 

As In hef. lO, the disk Is assumed to be rigid 
o."d the rotor Is modeled as an Infinite, two- 
dimensional cascade of airfoils capable of plunging 
and pitching motions. The geometry of a tuned cas- 
cade is shown In Fig. 1. As Illustrated In Fig. 2, 
each airfoil Is suspended by bending and torsional 
springs, Kn and Kq respectively. The 
blade Is assumed to be rigid In the chordwise direc- 
tion, and, consequently, this motion Is neglected. 

The coupling between bending and torsion due to 
effects such as pretwist, shrouds, and rotation of 
the rotor is modeled thruugh the offset distance (b 
x„ ) between me center of gravity and elastic 
axis positions. The centrifugal stiffening effects 
due to rotation are Included In the bending and tor- 
sional spring constants. The properties of the 
blades are represented by their respective values at 
75 percent span. 

In the case of a cascade with identical blades 
(tuned). It is usually assumed that the motion of the 
blades Is simple harmonic and that they vibrate In 
each mode of the cascade with a constant amplitude 
and with a constant phase angle between adjacent 
blades. The phase angle Is restricted to N 
discrete valies, 6r - 2«r/M, where r « 0, 1, 

2 ... N-1, and thus the total number of allowed 
Interblade phase angle modes of the cascade is also 
N. This restriction is known as Lane'sl^ assump- 
tion. In the case of a mistuned cascade considered 
herein, the general motion of a blade can be ex- 
pressed as a sum of Its motions In all possible 
Interblade phase angle modes. Then, the motion of 
the sth blade Is 



For a cascade with N mistuned blades Eq. (1) can be 
generalized as 

{X} e^"‘ - [E]{Y} e’“‘ (2) 


where 
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B. Aerodynamic Model 

As mentioned earlier, the unsteady aerodynamic 
loads were calculated by using Smith's^* theory In 
subsonic flow, and Adamcryk and Goldstein ' s*^ the- 
ory In supersonic flow with a subsonic leading edge. 
In these theories, the effect of airfoil thickness, 
camber, and steady state angle of attack are neg- 
lected, and the flow is assumed to be Isentropic and 
Irrotatlonal. In the calculation of subsonic aero- 
dynamic loads some numerical convergence problems 
have been encountered when Bp ■ 0 or Since 

the Br ■ 0 n*' mode Is of little Interest 1)i 
practice and since a numerical study of this problem 
Is beyond the scope of the present paper, this 
numerical problem nas been resolved by setting 
Br • O.OOl whenever Its value should be 0 or 2«. 

The effect of wakes shed from upstream periodic 
obstructions are Included. The wakes considered are 
limited to sinusoidal distortions represented by 
vorticity perturbations so that they are convected 
downstream at the flow velocity V. It should b'? 
noted that Ref. 12 does not Include the effect rf 
wakes from upstream periodic obstructions. However, 
the Inclusion of this effect In the formulation and 
the details of the computer program In the supersonic 
flow have been obtained from Adamczyk.l^ 

C. Equations of Motion 

The equations of motion for the cascade were 
developed In Ref. 10. For completeness, those equa- 
tions are summarized below without the derivation de- 
tails. The equation' of motion for the sth blade are 
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[A] 


° [A,] 


Ws “ fnj/tpb^ 


{AO}. [{ADpUAOj} ... {AD^.^}]^ 



D, Possible Regimes of Flutter 

It is obvToiis that in incompres*' Ible flow there 
are no acoustic (pressure) waves sMice the governing 
equation for the perturbation pressure or velocity 
potential is a Laplace equation. However, in com- 
pressible flow acoustic i/aves are possible since the 
equation for ttie perturbation pressure is a 'wave' 
equation. The possibi’ity of a wave type solution 
depends on the flow parameters. Different regimes of 
flutter 'or subsonic flow were categorized in Ref. 15 
based on whether an acoustic disturbance decays or 
propagates. The marginal condition between these two 
cases is called the 'acoustic resonance' or 'cut- 
off'i*^ condition. For more details see discussions 
in Refs. 11 and 15. A discussion was presented in 
Ref. 17 on different regimes of flutter for super- 
sonic flow. To facilitate later discussion of re- 
sults, a brief description of the possible regimes of 
flutter in both subsonic and supersonic flow follows. 


Subsonic flow . The decaying or propagating 
nature of an acoustic disturbance car. be described by 
the equation 




sin 



(7) 


gate. This marginal condition is called 'acoustic 
resonance* condition. If flutter occurs very close 
to acoustic resonance condition, it is cal’ed 
'acoustic resonance flutter.' The dutinction be- 
tween these regimes of flutter is not alw/.ys clear 
cut because they may merge, depending on the system 
parameters. Also, it should be remarked that this 
categorization of different regimes of flutter is 
mainly for descriptive purposes. The value of the 
integer n in Eq. (7) for the given values of H, k, 
s/c, and t should be determined such that the roots 
of the quadratic equat'on (obtained with equality 
sign in Eq. (7)), 6^1 and 6 f 2 , satisfy the con- 
dition, 0 < Brl 2 <_ 2*. Alternatively, sub- 
critical and supercritical regimes of flutter can 
also be categorized based on the values for Srl 
and 8r2* Flutter is called supercritical it the 
flutter interblade phase angle lies between Brl 
and fjf2 and is called subcritical otherwise. 

Su personic flow . The propagating nature of an 
acoust ir~disturbance in supersonic flow is also 
governed by Eq. (7). The acoustic resonance condi- 
tion obtained from Eq. (7) is the same as the 
Eq. (A 8 ) of Ref. 18. Also, the same resonance condi- 
tion can be derived from Eq. (1) of Ref. 17, In con- 
trast to the subsonic case, in supersonic flow when 
the left hand side of Eq. (7) is greater than zero, 
the acoustic disturbances propagate through the en- 
tire flow but remain bounded in the far field. If 
the flutter occurs when this condition is met, it is 
called subcritical flutter. If the right side of 
Eq. (7) is less than zero, the disturbances attenuate 
in the far field. The flutter, if occurring in this 
region, is called superciitical . 

III. SOLUTION 

The aeroelastic stability of the cascade is 
determined by the eigenvalues, y's, of the matrix [P] 
in Eq. ( 6 ). The relation between the frequency u 
and Y IS 


iu/u^ “ i/v^ = w * iv 


( 8 ) 


Flutter occurs when IT > 0. For the given values of 
the number of blades, and hence the allowable 
Br’s. the gap to chord ratio, the stagger angle, 
the elastic axis position, and the structural para- 
meters, the eigenvalues of the matrix [P] are cal- 
culated for a range of values of k. Denoting the 
values of k ancl ~ at which li = 0 as kp and 
vp , respectively, the nondimensional flutter speed 
can be written as 


Vp/buo ® Vp/kp 


(9) 


which is derived from Eq. (11) of Ref. 11. If the 
left hand side of Eq. (7) is greater than zero, the 
disturbances decay exponentially with distance from 
the cascade and wave type solution for the equation 
for perturbation pressure is not possible. If flut- 
ter occurs when this condition is met, it is called 
'subscritical flutter.' If the left hand side of 
Eq. (7) is less than zero, a wave type solution for 
the equation for perturbation pressure is possible 
and at least one pair of acoustic waves can propa- 
gate. If flutter occurs when this condition is met, 
it is called 'supercritical flutter.' If the left 
hand side of Eq. (7) is zero, a pair of acoustic 
waves is just on the verge of being able to propa- 


The aeroelastic response or the blades induced by 
wakes is calculated from Eq. (5) and is 

(Y) . - [[P] - [ijY]-’[Ej-iLG3L:3{AD} (10) 

IV. RESULTS AND DISCUSSION 

A. C o mputer Program and Verification 

To calculate the Hutter boundaries and the 
blade aeroelastic response of an arbitrarily mistuned 
cascade, a digital computer program was written based 
on the present formulation. In this program, it is 
possible to consider any type of mistuning such as 
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bUde to blade variations of the olade uncoupled 
bending and torsional frequencies, damping ratios, 
mass ratios, elastic axis and center of gravity posi- 
tions, and so on. This program Is operational on the 
NASAoLewls Research Center IBM 370/3033. Both the 
tuned and mistuned uncoupled bending and uncoupl’'d 
torsion cases, in addition to the tuned coupled 
bend Ing-tors Ion case, can be treated as special caies 
of this program. This program was checked for the 
following special cases: 

1. The subsonic aerodynamic coefficients were 
checked by comparison of the present results to the 
published results In Refs. 6 and 11. The supersonic 
aerodynamic coefficients were checked uy comparison 
of the present results to the published results In 
Refs. 12 and 18. 

2. To check the correctnesi of the program In 
predicing uncoupled torsional flutttfi of a tuned cas- 
cade In sub-:or 1c flow, a l?-bladsd rotor described In 
Ref. 1, was analysed, A comparison was made of the 
present results obtained as a special case uf the 
program to those presented In Ref. 16 for the 
l2-b1aded rotdr. 

3. The correctness of the program In predicting 
the coupled bend Ing-tors Ion flutter speed of a tuned 
cascade In supersonic flow was checked by comparing 
the present results for a few selected cases to the 
corresponding ones In Ref. 19. 

4. The correctness of the program in calculating 
coupled bending-torsion flutter and response of a 
mistuned cascade for M • 0 was checked by comparing 
the present results with the corresponding ones in 
Ref. 10, 

In all of these cases good agreement was 
obtained. 

B. Aetoelastlc Stabi'lty 

BotH the aeroelastic stability and response 
analyses presented in this paper are for NASA Test 
Rotor 1? (shown in Fig. 3). This special test rotor 
is similar to a forward stage of an advanced axial 
flow compressor. The required parameters of this 
rotor are listed in Table 1. It should be mentioned 
that this rotor was analysed for incompressible flow 
in Ref. 10. Itie blade bending-to-torsion freouency 
ratio for the tuned rotor is 0.367, and the elastic 
axis and c.g. position are at 60 percent chord. As a 
result, the coupling between bending and torsion Is 
very weak and the flutter mode Is dominated by tor- 
sional motion. Hence the results for the pre- 
dominantly bending modes for some cases will not be 
presented. However, to conduct parametric studies 
the bending-to-torsion frequency ratio and elastic 
axis position are varied. For this case the results 
for the predominantly bending mode will also be pre- 
sented. Furthermore, no attempt is made to match the 
flutter Mach number with the assumed Mach number in 
the parametric study results presented. 

As discussed earlier, three different regimes of 
flutter were found in both subsonic and supersonic 
flows. It is useful to illustrate them separately. 
Figures 4 to 6, show the variation of the real part 
of the eigenvalue, “ (which is a measure of stability 
and is defined in Eq. (8)), with the reduced fre- 
quency, k, at different Mach numbers for a tuned cas- 
cade. The flutter mode in these figures is dominated 
by torsional motion of the blades. In Fig. 4 for 
M « 0.86, the regions of subcritical and acoustic 
resonance flutter and the acoustic resonance reduced 
frequency, k^^. are shown for 8r - 83.57’, which 
is the critical interblade phase angle for sub- 
critical flutter. Also Indicated are the acoustic 
resonance flutter regions for a few other 6r's. 


No Supercritical flutter was found when M £ 0.B6. 

In Fig. 5 for M ■ 0.9, the regions of supercritical, 
acoustic resonance and subcritical flutter, and 
kjr are shown for - 61.45’ which Is the 
critical Interblade phase angle for supercritical 
flutter. The Intervals In which acoustic resonance 
flutter occurs for a few other values of 
shown. For a fixed a,, the subcritical flutter 
reduced frequency kp Is less inan the acoustic 
resonance flutter kp, which is less than the 
supercritical flutter kp. The appropriate value 
of n in determining the regimes of flutter In 
Figs. 4 and 5 is zero. By comparing Figs. 4 and 5 It 
can be x len that as N increases toward unity, tne 
regions of acoustic resonance occur at decreasing 
values of k. The frequency separation of these re- 
gions al'o decrease. A similar behavior was observed 
In Ref. 16 In which only uncoupled torsional motion 
Is considered. In Fig. 6 for M - 1.1, the regions 
of supercritical, acoustic resonance, and subcritical 
flutter are shown. For each there Is a range 
(a band) in which the supercritical flutter occurs. 
The width of the band depends on Bri And the bands 
for some B^'s overlap. The appropriate value of 
n Is negative one. No supercrit .cal flutter is 
found when the Mach number is greater than 1.15. 

Thus It Is hypothi'sized that supercritical flutter 
only occurs when \he Mach number is near unity 
(e.g., 0.85 < M < 1.16). 

A comparison of the system eigenvalues of both 
tuned and mistuned cascades is useful to understand 
the effects of mistuning on flutter and response. 
Figure 7 provides such a comparison for M • 0. As 
mentioned earlier, for xhe parameters considered in 
this figure the flutter mode is dominated by tor- 
sion. As a result, the eigenvalues corresponding to 
predominantly bending motion are of little interest 
and are not shown. The type of mistuning considered 
is the one in which the odd and even numbered blades 
have different uncoupled torsional frequencies. This 
is known as alternate blade mistuning. For example, 
in the case of 1 percent mistuning, the frequency 
ratio, u>a /uiQ, is 1,005 for all the even blades 
and is 0.996 for all the odd blades. The reference 
• 'equency wq is equal to the arithmetic mean of 
the uncoupled torsional frequencies of all blades. 
Because of the symmetry of this type of mistuning 
the 6r mode couples with the (Br * *) mode 
only. It should be pointed out that the value of ir • 
reduced frequency in this figure is chosen such that 
the cascade is neutrally stable for 2 percent alter- 
nate mistuning. Two interesting observations follow 
from Fig. 7. First, for M - 0 the 2 percent mis- 
tuning significantly affected the system eigenvalues 
and stabilized the unstable tuned cascade. As the 
level of mistuning is increased, the spread of the 
real part of the eigenvalue is decreased. This ob- 
servation implies that the effective damping of some 
modes is increased while that of other modes is de- 
creased with an ncrease in the level of mistuning. 
This behavior will be recalled in the discussion of 
forced response which will be addressed later. 

Second, as mistuning increases the modes separate 
into high and low frequency groups. The degree of 
separation increases with the level of mistuning. 
These findings art in agreement with the similar ones 
in Ref, 10. 

To illustrate the effect of Mach number in the 
presence of mistuning, the analyses presented in 
Fig. 7 are repeated in Figs. 8 to 10 at Mach numbers 
0.5, 1.15, and 1.4, respectively. Comparing the 
ranges of the real parts of the eigenvalues in these 
figures, several interesting conclusions follow. 
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First, this type of mistuning elMays has a stabiliz- 
ing effect for the Mach numbers considered. However, 
the stabilizing effect is stronger in subsonic flow 
than that in the supersonic flow. Second, in sub- 
sonic flow the range of the real parts of the eigen- 
values increases with increase in Mach number. But, 
ir the supersonic flow the opposite is true. When 
M a 1.4, the effective aerodynapnic damping in all the 
modes is relatively small. Third, the Mach number 
has a significant effect on the amount 'f mistuning 
required to separate the modes into high and low fre- 
quency groups. For example, Z percent mistuning 
separated the modes in supersonic flow but did not in 
subsonic flow. 

Figure 11 shows how the reduced frequency kp 
depends upon Mach number, mistuning, and structural 
uamping. Since the unsteady cascade aeiodynamic 
theories used herein are not valid in the transonic 
region, no attempt is made to extrapolate the curves 
in this region. Also, the acoustic resonance f'utter 
regimes are not shown. It is seen that as the Mach 
number increases from zero, the kp of the tuned 
cascade decreases until supercritical flutter 
occurs. This decrease in kp corresponds to 
either higher fluid velocities or lower blade stiff- 
ness. The effect of increasing M is, therefore, 
highly favorabW on coup’ed bemJing-torsion flutter. 
In contrast, when M increases from 0.86 to 0.9, its 
effect is unfavorable, and the flutter is of super- 
critical type. When M increa'.es from 1.1 to 1.1b, 
its effect is again favorable and the flutter is of 
supercritical type. This increase in kp corre- 
sponds to either lower fluid velocities or higher 
blade stiffness. Thus, the effect of increasing M 
on tbe kp of a tuned cascade depends on the range 
of M. Also, it is evident from Fig. 11, that the 
variation of kp wilti M in an alternately mis- 
luied tasiade depends on both the level of mistuning 
and 111 ' range of M. Ihe effect of structural dantp- 
iiK! 1 * I tuned cascade also depends on the range of 

> 1 1 iK t lira 1 ilaniping is highly effective when M 
!■. k'.s tian 0,6 and l»’ss effective when M is 
, I eater ttuin 0.6. I*, subsonic flow all levels of 

mistuning considered have much stronger effect on 
the kp than does O.i? percent structural damping. 
Whereas, in supersonic flow when M is between 1.16 
and 1.,^‘j, the 0.^’ percent damping is more effective 
than i’ percent mistuning and is less effective than 
5 percent mistuning. Wlien M is between l.<?b and 
!.b, the O.Z percent damping is more effective than 
b percent iiusluning and less effective than ,?0 per- 
cent mistuning. When M > 1.3, t.'ie Q,'i percent damp- 
ing is more effective than ?0 percent mistuning. 

It should be mentioned that predicted flutter 
Mach numbers for NASA Test Rotor 12 can be obtained 
by constructing an 'operating line' in Fig. 11. This 
operating line was constructed by fixing the speed of 
sound and ttiu blade uncoupled torsional frequency. 

The intersection of this line with the neutral sta- 
bility boundary curve determines the flutter Mach 
number. Although this line is not shown in fig. 11, 
the resulting flutter Mach number is 1.37b for the 
tuned undamped case. Of course, this value changes 
if one includes structural damping and blade mistun- 
ing. The flutter Mach numbers are 1.396 for 2 per- 
cent alternating mistuning, 1.415 for 5 percent mis- 
tuning, 1.440 for 20 percent mistuning, and 1.469 for 
0.2 percent damped tuned case. All of these flutter 
Mach numbers are above the design value. 

It is apparent from Fig. 11 that mistuning has a 
greater influence in subsonic than In stipersonic 
flow. To determine if this trend shown by Rotor 12 
is also similar for a fan stage, a cascade represen- 


tative of an advanced unshrouded fan stage (aspect 
ratio 3.34) was analyzed. It was found that mistun- 
ing had a strong effect on stability at a tip rela- 
tive Mach number of 1.32. This effect was similar to 
that shown by Rotor 12 at low subsonic Mach numbers. 
Additionally, it appears that the stabilizing effect 
is sufficient to utilize mistuning as a passive '.on- 
trol for unshrouded fan designs in supersonic flow. 

A detailed description of this analysis will be pre- 
sented in a future publication. 

The effect of the bend ing-to-tors ion frequency 
ratio, wb./wo< 4nd elastic axis position on 
flutter s^ed, Vp, is of general interest. To 
explore these effects, the parameters of the rotor 
are varied. For example. Fig. 12 shows the depend- 
ence of Vp of a tuned cascade (wo - i«o 
for all s) upon whc/"o The elastic axis 
at 25 percent chord with M as a parameter. When 0 
< wb /un < 0,35, the Mach number has favorable 
effect in“both subsonic and supersonic flows; when 
0.3' < wbr/wo < 1.5, M has unfavorable effect 
in supersonic flow and favorable effect in subsonic 
flow; and when I uq > 1.5, M hss unfavor- 
abl'? effect in both subsonic and .upersonic flows. 

A very limited study was conducted to determine 
the effects of alternating structural damping mistun- 
ing on aeroelastic stability. The odd numbered 
blades were left undamped and 0.4 percent structural 
damping was added to the even numbered blades. For 
this case with M » 0.5, the reduced flutter fre- 
quency kp is 0.517. As can be seen front Fig. 11, 
this value is very close to tl o kp - Q.511 for the 
case where 0.2 percent daiiipinc vas added to all the 
blades equally. Also from Fig U, the value of 
kp for the undamped tuned cascade is 0.625. 
Comparison of these three values of kp shows that 
the alternating structural damping does not have any 
additional benefit over adding damping to all blades 
equally. Comparison of the eigenvalues, which are 
not shown due to space limitations, further shows the 
alternating damping mistuning does not result in the 
significant mode coupling as was noticed for alter- 
nating frequency mistuning. 

C. Aeroelastic Response 

In tlie present formulation, it is possible to 
consider an excitation cons'sting of all harmonics of 
rotational speed of the rotor which range up to 
N-1. The coefficients, ly,hr *nd Iwar. I'’ 
to. (4) represent the forcing functions in the bend- 
ing and torsional equations, respectively, in the rth 
mode. To understand the nature of the response, 
excitation in only one harmonic at a time will be 
considered. This results in no loss of generality 
since the principle of superposition holds. If the 
r - R harmonic is considered, then the column 

matrices {ADq}, {ADj} {ADn_i) are zero 

except ADp in Eq. (6). According to the travel- 
ing wave representation in Eq. (1), this corresponds 
to N-R symmetrically space obstructions located 
upstream from the b’xdes and the circumferential wake 
distribution is perfectly sinusoidal. For practical 
applications, the forcing frequency is thus equal to 
(N-R) times the rotational speed. 

The aeroelastic response is presented for two 
values of R, 11 (45 obstructions) and 39 (17 ob- 
structions), at a fixed reduced frequency which is 
chosen such that the cascade is aeroelastical ly 
stable in all modes. These values for R were 
picked because the aerodynamic damping of the tuned 
system in the r . 11 mode is relatively low, whereas 
th-jt in the r - 39 mode is relatively high. The 
forcing frequency range investigated is limited to a 
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small ranae around th« uncoupled torsional fre- 
quency. If the blades are tuned, the response will 
be entirely in the r - R mode, and all the blades 
have equal amplitudes. The amplitude of response Is 
a function of w/wq. If the blades are now arbi- 
trarily mistuned, there will be response In all the 
inodes (enumerated by r) and the response of the sth 
blade can be obtained from Eq. (1). Figure 13(a) 
for R • 11 and Fig. 13(b) for R ■ 39 show the 
variation of aj/aj for both tuned and the 
two percent alternaling blade mistuned cascades In 
subsonic flow. The quantity as,1d The ampli- 
tude of the torsional motion of each blade at reso- 
nance In the tuned rotor and it depends upon R. 
Figures 14(a) for R •• 11 and 14(b) for R « 39 are 
repetitions of Figs. 13(a) and (b), respectively, in 
supersonic flow. The bending amplitudes are not 
shown because they are very small In the range of the 
excitation frequency shown herein. For alternating 
mistuning, only the Sr (^r * ■) modes are 

coupled. Therefore, In all three cases the single 
resonance peak of tne tuned cascade Is replaced by 
twin resonance peaks- It Is seen from the figures 
that the effect of mistuning on forced response de- 
pends upon the engine order of the forcing function 
and upon the Mach number. For example, In subsonic 
flow the mistuning has a beneficial effect on tor- 
sional response (Fig. 13(a)) for R > II, but has an 
adverse effect (Fig. 13(b)) for R - 39. This be- 
havior Is similar to that noticed in the incompressi- 
ble flow In Ref. 10, But, this behavior 1s In con- 
trast to the common belief thct mistuning always has 
an adverse effect on forced response. Therefore, 
this observation provides an added incentive for 
pursuing the use of mistuning as a passive control. 
The decrease in the ma*iniui.’ amplitude of any blade 
with mistuning for R « 11 1, approximately 51 per- 
cent (fig. 13(a)) and the increase in the maximum 
amplitude for any blade with mistuning for R » 39 Is 
approximately 60 percent (Fig. 1.1(b)). In supersonic 
flow for M - 1.1, the alternating mistuning has an 
unfavorable effect for both the values of R - 11 and 
39. One of the reasons for this behavior is that the 
tuned modes R » 11 and R r 39 do not iiave relative- 
ly low and high aerodynamic dampings as in the sub- 
sonic flow. 

V. CONCLUSIONS 

The analyses of the effects of mistuning on tne 
coupled bending-torsion flutter and response for In- 
compressible flow which were developed in Ref. 10 
have been extended into the subsonic and supersonic 
flows. The following conclusions are reached based 
on the limited parametric studies presented In this 
paper. 

1. All three regimes of flutter (sohcritical , 
acoustic resonance, and supercritical) were found to 
exist in both subsonic and supersonic flow for cer- 
tain parametric combinations. However, supercritical 
flutter was found only when the Mach number Is near 
unity for the rotor considered herein. 

2 . It was found that mistuning can have a sig- 
nificant effect on flutter speed depending on the 
cascade parameters. For Rotor 12 the potential was 
found to be greater in subsonic than In supersonic 
flow. For a supersonic fan the effect appear, to be 
sufficient to utilize mistuning as a passive flutter 
control . 

3. The use of alternating structural damping was 
not found to have any additional benef t over adding 
damping to all blades equally. 


4. The addition of a small amount of structural 
damping was found to have a much stronger effect for 
low subsonic Mach numbers than vor high subsonic or 
for all supersonic Mach numbers. 

‘ 5. As reported by the authors for Incompressible 

fUw, the use of uncoupled torsional flutter analysis 
to deduce the effect of elastic axis position was 
found to be unreliable. Coupling between bending and 
Version, structural damping, and iristuning can change 
the results significantly In both subsonic and super- 
sonic flows. 

6. In general an Increase In Mach number was found 
to have a favorable effect on subcritical flutter In 
subsonic flow and have an unfavorable effect on 
supercritical flutter In supersonic flow. Its effect 
also depends on the bending-to-torslon frequency ratio 
when the elastic axis Is off mid-chord. 

7. Mistuning may have either a beneficial or an 
adverse effect on forced response, depending upon the 
engine order of the excitation and upon the Mach 
number. 
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Figure 1. - Geometry of a tuned cascade in rthmode (note that the variable 
s represents the gap between the chords as well as the blade number 
indexi. 
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Figure 2. • Airfoil restrained from bending and torsional motion (si- 
nusoidal wakes not shown). 
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Figur* 5. - SubcrHIcjl. suptrcrltical. and acoustic rtsonanca fluttar 
predictions ot a tunad cascada: M • a9. a • ^ x„ • a wn '‘^o, * 

a»7. 
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Figure 6. • Subcritlcal. acoustic resonance, and : jpercritkal flut- 
ter predictions of a tuned cascade: M*L1, a*0, *(i 
• Ol 357, undamped. 




I 




o<in 


& 8 » S ^ 

— • • • 

^ *3mVAN39l3 30 IMVd AdVNIOVWI 



I 

I 


c W 

S o 
C xT 

•» ^ 

SIS 

8a 

i- 

•R cf 

i> 

If 

■53 !Q 
■8 -» 



i 






TORSIONAL AMPLITUDE RATIOS. TORSIONAL AMPLITUDE RATIOS, 

id id 


TUNED 



Figure 13(al. - Effect of blade alternating mistuning on coupled bend- 
ing torsion response: M ■ 0. 5, R - 11, Cf, • Cq • OL 002, 

- a 357. * * 



Figure 13(b). Effect of blade alternating mistuning on coupled bend- 
ing torsion response: M • 0. 5. R •> 39, ^ 



Figure 14(a). - Effect of blade alternating mistuning on coupled bend- 
ing torsion response-. M • 1. 1. R - 11, Cu • 0.002, 

Whj'-’o ■ a ^ 
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Figure 14(b). - Effect of blade alternating mistuning on coupled bend- 
ing torsion response: M • 1.1, R • 39, *Cc, ■ 0.002, 
too -a 357. ' ^ 



